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RAT J0455+1305: A rare hybrid pulsating subdwarf B star 
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ABSTRACT 

We present results on the second-faintest pulsating subdwarf B (sdB) star known, 
RAT J0455+ 1305, derived from photometric data obtained in 2009. It shows both short and 
long periods oscillations, theoretically assigned as pressure and gravity modes. We identify 
six short-period frequencies (with one being a combination) and six long-period frequencies. 
This star is the fourth hybrid sdB star discovered so far which makes it of special interest as 
each type of mode probes a different part of the star. This star is similar to the sdB hybrid pul- 
sator Balloon 090100001 in that it exhibits short-period mode groupings, which can be used 
to identify pulsation parameters and constrain theoretical models. 
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1 INTRODUCTION 

Hot subdwarf (sdB) stars are horizontal branch stars with 
masses near to 0.5 Mq and very thin (in mass) hydrogen 
envelopes. Average effective temperatures and surface gravi- 
ties are about 30000 K and \ogg ~5.5, respectively. Al- 
though it is clear that sdB stars will eventually enter the 
white dwarf cooling track witho ut reaching the AGB phase 
jPorman. Roo d & 0'Conneiil ll993'), their formation as sdB stars 
is under debate. There are sev eral mechanisms that involve single- 
star or binary evolution (e.g. D'Cruz. Dorman. Rood et al.|[l993 : 



terest are these rare stars showing both (p— ) and (y— )modes 
since models of these stars can constrain their cores and outer 
regions simultaneously. The first h ybrid sdBV star was found by 
ISchuh. Heber. Dreizler et all (|2006|) and three more have been dis- 
covered using ground-based data, with RAT J0455+1305 being the 
faintest. The Kepler satellite has also recently discovered several 
candidate hybrid sdBV star s (Reed. Kawal er. 0stensen et al.ll2010l: 



Kawaler. Reed, Quint et alJl2OlO, ; ,0stensen. Green. Bloemen et aU 
2010bh from its survey phase with follow-up data to be obtained 



iHan. Podsiadlowski. Maxted et alj|2002l. |2003|) . Detection of pul- 
sations in hot sdB stars opened a way to study their interi- 
ors and evolution prior to the horizontal branch. First, short pe- 
riod os cillations were found bv lKilkennv. Koen. O'Donoghue et afl 
( Il997h in EC14026 (now officially named V391 Hya). Ac- 
cording to theoretical models these pulsations are attributed 
to pressure (p— )modes and are driven in the outer part of 
the stars ICharpinet. Fontaine. Brassard et al.lll997l). S everal years 
after their discovery. Green. Fontaine, Reed et all (l2003h an- 
nounced another kind of variability in sdB stars. The longer- 
period variations were also identified as stellar oscillations. In 
this case however, they are attributed to gravity (g— )modes 
and originate deeper within these star s than the (p— )modes 
jFontaine, Brassard, Charpinet et al.ll2003l) . Both types share some 
overlap in the H-R diagrsim and so the same stellar models are 
appropriate for studying both kinds of pulsation. Of special in- 
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over the next several years. 

The brightest known pulsating subdwarf B (sdBV) star 
has a V magnitude of 11.8, t he typical brightness is about 
14.5 a nd the faintest are 17 dOstensen. Oreiro. Solheim et al.l 
l2010ah . Methods for constrai ning the pulsation modes include 
multicolor photometry dTremblav, Fontai ne. Brassard et al.l 
2006) . low - resolution spectroscopy ([Telting & Oste nsenI 
2004 I2OO6I ; iReed, Eggen, Harms et all \200^. a co mbina- 
tion of these two l lBaran. Pigulski & 0'Tool3 l2008h . and 
time-series observations o f line profile variations using high- 
resolution spectrosocpy (Telting. Geier, 0stensen et"aLl I2OO8I : 
pTeltin g, 0stensen, Oreiro et al. 2010). This last method requires 
large telescopes, where it is quite difficult to get time. To date, 
time-series spectroscopy has been obtained for only a few objects, 
all brighter than 14 mag. Because of the difficulty of obtaining 
data spanning relatively long periods of time, the faintest stars are 
usually limited to discovery observations taken in one or no filter 
at all (so-called white light). 

Single-color data, which are typically obtained to increase 
temporal resolution, are usually not useful for mode identification, 
which is necessary for constraining the models. To identify modes 
using broadband single color data one needs to assume trial values 
for free parameters in models and search for the best fit to observed 
frequencies. The large number of free parameters does not lead to 
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Table 1. Log of observations 



Date 


hours 


exposure 


filter 


site 


10 Nov 2009 


7.2 


30s 


none 


San Pedro Martir 


1 1 Nov 2009 


7.4 


30s 


none 


San Pedro Martir 


14 Nov 2009 


7.7 


30s 


none 


San Pedro Martir 


24 Dec 2009 


1.9 


30s 


BG40 


KPNO 


25 Dec 2009 


8.6 


30s 


BG40 


KPNO 


26 Dec 2009 


8.6 


30s 


BG40 


KPNO 


27 Dec 2009 


8.4 


30s 


BG40 


KPNO 



strong constraints on stellar models. An exception is when the fre- 
quency spectrum contains specific features, such as frequency mul- 
tiplets, mode groupings, or small and large spacings, which con- 
strain mode identifications and therefore reduce the number of free 
parameters. 

Following this approach, and encouraged by the results ob- 
tained for other sdBV stars, we decided to obtain photometric 
observations of one of the faintest sdBV stars discovered so far, 
RATJ0455+1305 (hereafter RAT0455). Our results clearly show 
that small and/or moderate telescopes are useful for the asteroseis- 
mic study of faint sdBV stars. 



2 RAT0455 

RAT0455 was discovere d during the RApid Temporal Survey 
( iRamsav & Hakalall2005l) . The aim of this survey was to search 
for variability from a few minutes up to several hours, for stars 
down to 22 mag to find interacting ultra-compact binary systems 
with orbital periods less than 70 min. The design of the project is 
also useful for detecting sdBV stars, which can serendipitously oc- 
cur in the field of view. To date, as a by-product of this survey, 
one sdBV star has been discovered. It was assigned the designation 
RAT J0455+ 1305 in accordance with the survey convention. Dis- 
covery observations along with Fourier analysis and a classification 
spectr um were published by iRamsav. Napiwotzki. Hakala et al.l 
( l2006h . RAT0455 is 17.2 mag in the the V filter and the ampli- 
tude spectrum calculated from the discovery data revealed only one 
frequency in the (p— )mode region with a relatively high ampli- 
tude. Placing RAT0455 on the logTcff and log g diagram (Fig[T}, 
we can see that it is s imilar to other hybrid stars, and indeed 
iBaran & Fox-Machadd l l2009l) detected peaks at low frequencies 
which were assigned as (g— )modes. 



3 OBSERVATIONS 
3.1 First data 

Our first effort to obtain follow-up time-series photometry of 
RAT0455 was in January 2009 using the 1.5 m telescope at San 
Pedro Martir Observatory in Baja, Mexico. To maximize sig- 
nal, the observations were conducted without a filter. Photomet- 
ric data were obtained during three of five nights. Unluckily, non- 
photometric conditions deteriorated the data quality. Those data 
confirmed the dominant periodicity at the frequency and ampli- 
tude reported in the discovery paper, and revealed three new small- 
amplitude peaks in the long period domain. The latter are usually 
identified as (g— )modes in sdB stars. Detection of both short and 
long period modes confirmed that RAT0455 is indeed a hybrid star 
and increased the number of known hybrid sdBV stars to 4. Details 



Table 2. Results of the pre-whitening process for KPNO data. The phases 
are given at mean epoch 2455145.89952. The numbers in parentheses are 
the errors of the last digits and the residual noise level is 0.26 mma. 



mode 


Frequency [mHz] 


Amplitude [mma] 


phase [rad] 


fA 


0.22930(21) 


2.2(2) 


1.53(10) 


fB 


0.24910(36) 


1.3(2) 


5.99(17) 


fc 


0.27837(22) 


2.1(2) 


3.44(10) 


fo 


0.30793(25) 


1.7(2) 


1.52(12) 


fE 


0.31398(21) 


2.1(2) 


0.92(10) 


fF 


0.48251(37) 


1.1(2) 


2.68(19) 


fl 


2.68044(2) 


18.8(2) 


0.90(1) 


h 


2.74899(10) 


3.9(2) 


4.25(5) 


h 


2.67673(9) 


5.3(2) 


5.45(4) 


f4 


3.59305(33) 


1.2(2) 


5.47(17) 


fs 


4.40179(36) 


1.1(2) 


0.61(18) 


fl+2 


5.42943 


0.8(2) 


5.50(24) 




log Teff 

Figure 1. Location of hybrid stars in the logT^fj - logg plane. Previ- 
ously known hybrid stars are marked with colors: Balloon 090 100001 -blue, 
DWLyn-red, V391 Peg-green. The location of RAT0455 is marked with 
'RAT'. 
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Figure 2. Photometry obtained during late 2009 at two sites: San Pedro 
Martir (left panels) and KPNO (right panels). Closer views are in the lower 
panels. 
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Figure 3. Window function of data taken at KPNO. Large diurnal aliasing 
indicates the single-site nature of tlie data. 



Table 3. Frequency ratios of consecutive (p— )mode regions for RAT0455 
and Bal09. The ratios were calculated using the highest amplitude frequen- 
cies in each group. The fourth group in RAT0455 was not detected. 



groups 


RAT0455 


Bal09 


1/2 


0.746 


0.745 


2/3 


0.816 


0.813 


3/4 




0.850 




Frequency [mHi] 

Figure 4. Amplitude spectra for the original data and after some prewhiten- 
ing steps for the KPNO data. Each panel indicates the number of frequen- 
cies removed from the original spectrum. 




2.65 2.7 
Frequency [mHz] 

Figure 5. Amplitude spectrum ai'ound the dominant peak. Location of two 
prewhitened peaks ai'e indicated by solid lines. The dashed line indicates the 
position of the remaining peak and is added to indicate a potential triplet. 



on observations, analysis and con clusions can be found in the paper 
bv lBaran & Fox-Machadol ( l2009h . 



3.2 Late 2009 data 

One of the conclusions of our previous paper was that more data 
were needed to detect more, lower amplitude frequencies (if any 
existed) and to search for other features, such as evenly spaced pe- 
riods and frequencies or multiplets. Such features may help to con- 
strain pulsation parameters, in particular the three quantities that 
describe the pulsation geometry and are usually denoted as n, I 
and m. In late 2009, we were allocated time on two telescopes: the 
1.5 m at San Pedro Martir and the 2.1 m at KPNO. At the former 
we had five nights while at the latter we planned 7 nights. Again 
the weather was not cooperative and the overall duty cycle was 
only about 60% of the available time. At San Pedro Martir we used 
a Marconi CCD with a square 3 by 3arcmin field of view. The 
exposure time was 30 sec with 7 sec read out giving a cadence of 
37 sec. As in January we did not use any filter. At KPNO we used 
an Apogee CCD also with a 3 by 3 arcmin field of view. We used 
exposure time of 30 sec, with a cadence of 32 sec, but here we used 
a BG40 filter to improve the contrast with the sky background. The 
details are given in Table[T] All data were calibrated for instrumen- 
tal effects and the brightness was extracted by means of PSF fitting 
and aperture photometry. Those data which were characterized by 
smaller scatter were taken for the final analysis. The data were cor- 
rected for color extinction with each night treated separately. 

Because of non-photometric conditions on 1 1 Nov 2009, those 
data were very noisy and not used in our Fourier analysis. They are 
not plotted in Figure |2l which shows data obtained at San Pedro 
Martir (left panels) and at KPNO (right panels). The upper panels 
for each site show all data, while the lower panels zoom in on a part 
of one night from each site. From this plot one can see the different 
quality of data obtained at the two different sites. Data from the 
1.5 m telescope are noisier, and the lack of filter results in a broader 
bandpass which reduced the pulsation amplitudes. 



4 AMPLITUDE SPECTRUM 



To derive amplitudes, frequencies and phases of the pulsations 
a Fourier analysis and pre- whitening process were employed. At 
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each step of pre-whitening, sinusoidal terms were improved by 
means of a non-linear least-squares method. Since the observations 
in San Pedro Martir were performed without any filter, and were 
separated in time from the KPNO run, data from different sites were 
analyzed separately. 

4.1 San Pedro Martir data 

We only obtained useful data on two nights from San Pedro Mar- 
tir and because of a flux discontinuity, part of 14 Nov data were 
removed. This results in a low duty cycle for these two nights of 
data, with an extra gap in the second night. Combining these two 
nights causes an aliasing pattern of 1/(4 days) = 2.89 /iHz from any 
real peak which complicates the pre-whitening process. Analysis 
of the combined data revealed a timing issue in the data from 14 
November. It was caused by incorrect times being written to the 
FITS header. Unfortunately, the reason for this failure remains un- 
known and so these data were not usable - we are therefore left 
with just one night's data. As such, these data did not detect any 
new frequencies. The only periodicity detected with S/N > 4 is the 
dominant peak at a frequency of 2.6782(7) mHz, an amplitude of 
13.1(4) mma and a phase of 4.19(3) rad at epoch 2455145.89952. 
This poor estimation of the frequency may be caused by bad cov- 
erage during our last night. 



4.2 KPNO data 

In December 2009, observations were obtained at the 2.1 m tele- 
scope at KPNO. Of seven allocated nights, data were collected dur- 
ing the first four. In total, more than 27 hrs of data were gathered. 
Because of the larger aperture at KPNO compared to San Pedro 
Martir, these data show a higher signal-to-noise, and even more im- 
portantly, no timing issues were present. Data from all four nights 
were combined and analyzed by means of Fourier analysis. Fig- 
ure [3] shows the window function calculated from these data. Fre- 
quencies, amplitudes and phases are presented in Table |2] Figure 
|4] shows residual amplitude spectra after a few steps in the pre- 
whitening process. The frequency resolution for these data is about 
3.5 /iHz. 

During these observation s, we detected more frequencies than 
iBaran & Fox-Machadd ( f2009l) did. Some were previously detected 
and are confirmed with these data while others are detected for the 
first time. Unfailingly, i t is easy to say that our domin ant peak fi 
is the same detected bv lBaran & Fox-Machadd ( l2009h and our f2 
is equivalent to theirs. For other frequencies, it is not as easy to 
compare the two investigations. Peaks f4 and fs from our anal- 
ysis were not detected before. Their locations are relatively reg- 
ular and similar to those found in Balloon 090100001 (Bal09) 
tearan. Pigulski. Koziel et alj2005h . In the (j-mode region, it is dif- 
ficult to compare peaks between seasons because the spectrum is 
very dense, so we do not know if two peaks with similar but not 
identical frequencies are the same (to within the errors) or if they 
are completely different peaks. Frequency proximity may suggest a 
similarity between fe and fs or fs and f^, where the former peak is 
detected in this analysis while the latter in lBaran & Fox-Machadd 
(12009'). This results in either 5 or 9 independent ((?— )modes de- 
tected so far. 

Also, we cannot investigate amplitud e stability as these data 
were not obtained in the same filter as iBaran & Fox-Machadd 
( l2009l) . so amplitudes cannot be directly compared. Nevertheless, 
they are for the dominant peak, i.e. 18.8 mma and 17.1 mma (in 



the previous season). Though the newest data were taken with the 
BG40 filter, we can conclude that the amplitude has decreased. Am- 
plitudes are decreased when no filter is used as it is averaged out 
through the wider bandpass. In total we detected 5 high-frequency 
peaks indicative of p-modes and 6 peaks in the low-frequency g- 
mode regime. We also detect one combination frequency of the 
dominant mode and f2. Its amplitude is very small and below the 
S/N > 4 criterion, but as a combination peak it is also included in 
the solution. 

What may be surprising is why, in the lowest panel in Figure 
|4l the residual amplitude spectrum still has a peak with relatively 
high amplitude which satisfies the above criterion. Its frequency 
falls near 2.6838 mHz and pre-whitening of this peak was unsuc- 
cessful. Including this peak grossly degenerated the fits to other 
peaks in the region. It should be noted that fi is separated from 
fg by approximately the formal frequency resolution, and this may 
cause unreliable fits via non-linear least-squares method. The same 
problem plagues the periodicity that we cannot pre-whiten. Its sep- 
aration from the dominant peak is almost the same as for fs. We 
suspect that with the current temporal coverage we cannot fully 
resolve all of the peaks around the dominant peak. As such, we 
emphasize that frequencies with small amplitudes in the vicinity of 
fi must be considered with caution. 



5 DISCUSSION - PRELIMINARY MODE 
IDENTIFICATION 

Using our detected frequencies, we would like to correlate them 
with the three pulsation parameters which describe pulsation ge- 
ometry, i.e. / and m, and the radial order n. Mode identification 
is essential for constructing theoretical models of stellar interiors. 
Experience from other pulsating sdB stars show that reliable mode 
identifications are not easy. 

These stars are faint, and that is why observations typically 
only use one filter. In such cases, identifying features like multiplets 
can help to deduce one or more pulsation parameters. Even nar- 
rowing the range of possible values can be an advantage to provide 
tighter constraints on theoretical models. Unfortunately, multiplets 
are not common among sdBV stars (Reed 2008). Usually peaks do 
not appear in a systematic way and then the only way to probe stel- 
lar models is to test each frequency against reasonable values of 
various free parameters. A solution is then attributed to the best fit 
between the observed frequencies (and global stellar parameters) 
and the model values jCharpinet, Fontai ne, Bra ssard et al..2005.) . 

In the case of RAT0455, a dozen frequencies have been de- 
tected. Luckily, in the amplitude spectrum there are features which 
can help in mod e identifications. They are similar to what was 
found in Bal09 jBaran. Pigulski. Koziel et al]|2005l) . First is the 
groupings of (p— )modes. These are not spread over a wide range 
of frequencies but instead they are located in specific regions only. 
In the case of Bal09 they were interpreted as consecutive overtones 
of the radial modes. This means the adjacent groups contain modes 
which differ in n by 1 . In RAT0455 such groupings are also present. 
However, since the noise level is higher, we could detect only one 
peak in each of two groups which satisfy the S/N > 4 rule. In the 
residual amplitude spectrum (lowest panel in Figure|4) one can see 
a signal excess near to fi and fs. It is likely that those groups con- 
tain more peaks but with amplitudes too small to be detected with 
our data. 

The frequency ratio between consecutive regions (calculated 
using the highest amplitude frequency in each group) as well as 
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those for Bal09 for comparison are presented in Table |3] Our re- 
sults are very similar to Bal09 which suggests that the interpretation 
for Bal09 may apply to RAT0455. Therefore, the dominant peak 
could be a radial fundamental mode, while the second and third 
groups of frequencies might contain consecutive radial overtones. 
This means that if fi is the radial fundamental (n = Z = 0) then f4 is 
the first overtone (n = 1, 1 = 0) and fs is the second overtone (n = 2, 
Z = 0). Of course, we could make a mistake assigning the dominant 
frequency to the fundamental mode instead of a higher overtone 
with radial order of n. Then fi and f 5 would have n = n+l and 
n = n+2, respectively. This includes an assumption that the highest 
amplitude frequency within each group is a radial (1 = 0) mode, as 
they suffer the least from surface cancellation effects. It does not 
have to be true, but mode identifi cations of other sdBV stars shows 
this to be the most likely case jBaran. Pigulski & O'Toolell2008l : 
iTelting & 0stensenll2OO4l) . 

Another feature even more useful for mode identifications is 
the existence of multiplets - three or more peaks equally spaced 
in frequency as a result of rotation (similar to the Zeeman effect 
in atomic structure). With multiplets, we can almost instantly as- 
sign I and m values to each component. In addition, a rotation 
period can be derived from the splittings. However Table |2] does 
not reveal any obvious multiplets. Figure 5 shows the region sur- 
rounding the dominant peak after removing fi, f2 and fs. It is ob- 
vious that there is another peak satisfying S/N > 4 rule, however 
its removal was not possible as it leads to spurious solutions (see 
I4.2t . It is possible evidence for a multiplet, but the insufficient fre- 
quency resolution does not allow us to resolve such closely spaced 
peaks. The two identified (and removed) peaks along with the one 
in the residual amplitude spectrum may form an almost equidis- 
tant triplet with a mean splitting of Af=3.5/xHz. If we assume 
that it is a rotationally split 1=1 mode then the rotation period is 
frot = (1 — Cn.i) I A/. Assuming that the Ledoux constant (Cn.i ) 
for Z = 1 is small (<0.05 lcharpinet. Fontaine. Brassard et aI.l20"ool) . 
we derive Prot to be about 3.3 days. Of course if these peaks are 
components of a triplet, then the dominant peak must have a non- 
zero value of I (i.e. I = 1). If so, this is contrary to the conclusion 
from the previous section that this peak may be radial. Additionally, 
it would mean that no radial modes are observed, despite their lack 
of any geometric cancellation. Another possibility is these three 
peaks do not form a triplet, but rather only two of them i.e. fs and 
the residual peak are parts of a triplet with a third member below 
our detection limit. The splitting would then be around 7 /^Hz and 
the corresponding rotation period close to 1.65 days. 



6 CONCLUSIONS 

We presented analysis of photometric data for one of the faintest 
sdBV stars, RAT0455. Despite its faintness, we obtained sufficient 
data to detect more, previously undetected frequencies and con- 
firm the hybrid nature of RAT0455. Similar to Bal09, we find non- 
uniform groupings of short-period p-modes and like Bal09, assign 
them to be the radial fundamental mode and two overtones. Such 
mode constraints apply directly to stellar models through which we 
learn about stellar interiors. A possible triplet is found, but higher 
frequency resolution is needed to confirm its existence. Our results 
do not, at this point, provide high confidence in our mode identi- 
fication , but this star looks very promising. Longer data coverage 
with a lower noise level in the amplitude spectrum would likely 
detect additional pulsation frequencies and answer these questions. 



That would make RAT0455 one of the very few sdBV stars with 
observationally constrained mode identifications. 



ACKNOWLEDGMENTS 

AB gratefully appreciates funding from Polish Ministry of Science 
and Higher Education under project No. 554/MOB/2009/0. JTG 
was supported by the Missouri Space Grant Consortium, funded by 
NASA. LFM acknowledges financial support from the UNAM via 
PAPIIT grant INI 14309. 



REFERENCES 

Baran A., Pigulski A., Koziel D. Ogloza W., Silvotti R., Zola S., 

2005, MNRAS, 360, 737 
Baran A., Pigulski A., O'Toole S., 2008, MNRAS, 385, 255 
Baran A., Fox-Machado L., 2009, Ap&SS, tmp, 285 
Chaipinet S., Fontaine G., Brassard P., Chayer P., Rogers F.J., 

Iglesias C.A., Dorman B., 1997, ApJ, 483, L123 
Charpinet S., Fontaine G., Brassard P., Dorman B., 2000, ApJ, 

131, 223 

Charpinet S., Fontaine G., Brassard P., Green E.M., Chayer P., 

2005, A&A, 437, 575 
Dorman B., Rood R.T., O'Connell R.W., 1993, ApJ, 419, 596 
D'Cruz N.L., Dorman B., Rood R.T., O'Connell R.W., 1996, ApJ, 

466, 359 

Fontaine G., Brassard P., Charpinet S., Green E.M., Chayer P., 

Billeres M., Randall S.K., 2003, ApJ, 597, 518 
Green E. M. et al., 2003, ApJ, 583, L31 

Han Z., Podsiadlowski Ph., Maxted P.F.L., Marsch T.R., Ivanova 
N., 2002, MNRAS, 336, 449 

Han Z., Podsiadlowski Ph., Maxted PEL., Marsh T.R., 2003, MN- 
RAS, 341, 669 

Kawaler S.D. et al., 2010, MNRAS, in press 

Kilkenny D., Koen C, O'Donoghue D., Stobie R. S., 1997, MN- 
RAS, 285, 640 

Ostensen R.H. et al., 2010, A&A, 513, 6 

Ostensen R.H. et al., 2010, MNRAS, tmp, L126 

Ramsay G., Hakala P, 2005, A&A, 360, 314 

Ramsay G., Napiwotzki R., Hakala P, Lehto H., 2006, A&A, 371, 
957 

Reed M.D., 2008 in U.Heber, S.Jeffery & R.Napiwotzki, eds. 
Hot Subdwarf and Related Objects ASP Conference Series, San 
Francisco, 392, 35 1 

Reed M.D. et al., 2009, A&A, 493, 175 

Reed M.D. et al., 2010, MNRAS, in press 

Schuh S., Heber U., Dreizler S., Heber U., O'Toole, S.J., Green 

E.M., Fontaine G., 2006, A&A, 445, 31 
Telting J. H., Ostensen R. H., 2004, A&A, 419, 685 
Telting J. H., Ostensen R. H., 2006, A&A, 450, 1 149 
Telting J. H., Geier S., Ostensen R. H. et al., 2008, A&A, 492, 815 
Telting J. H., Ostensen R. H., Oreiro R., Heber U., Vuckovic M., 

Randall S., Baran A., 2010, Ap&SS, tmp, 88 
Tremblay P.-E., Fontaine G., Brassard P., Bergeron P., Randall 

S.K., 2006, ApJS, 165, 551 



